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It is shown theoretically that the Sommerfeld coefficient j(B) in the mixed state contains useful 
information on the gap topology in a superconductor. We accurately evaluate 7(B) against magnetic 
induction B by using microscopic quasi-classical theory. The linear portion in 7 (ex B) relative to 
_B C 2 and its slope in low field faithfully corresponds to the gap anisotropy, so that we can distinguish 
the nodal gap and the anisotropic full gap, and estimate the magnitude of the gap anisotropy. 
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There has been much attention focused on exotic su- 
perconductors in recent years, because various new types 
of superconductors are successfully synthesized. For ex- 
ample, high T c cuprates and CeCoIn5 are thought that 
the superconducting state with <i-wave pairing is real- 
ized. In the case of UPt3 [lj and Si'2Ru04 |2|, the super- 
conducting state is thought to be a spin triplet pairing 
state. Only for a few superconductors, however, the re- 
alized pairing function has been conclusively determined 
as for the gap structure and its parity. Needless to say, it 
is crucial to know these when we understand its pairing 
mechanism. 

One of the standard methods to resolve the gap struc- 
ture is to use the temperature (T) dependence of phys- 
ical quantities, such as the specific heat C(T), the ther- 
mal conductivity k(T) and the spin lattice relaxation 
time Ti l (T). These comprise a unique set of charac- 
teristic T-dependences, allowing us to extract the gap 
topology; In the case of the nodal gap, where the gap 
in momentum space is zero along a line on the Fermi 
surface, C(T) cx T 2 , k(T) cx T 2 and Tf 1 (T) cx T 3 In 
contrast, these are all activation type for full gap case. It 
is often the case that these T-dependences are masked by 
other contributions such as impurities, phonon, nuclear 
spin, etc. Thus it is difficult to definitely conclude the 
gap anisotropy by the T-dependence of these quantities 
alone. We certainly need another methods to supplement 
this situation on top of existing ones. 

The magnetic field (B) in a superconductor adds a 
new dimension: The specific heat and thermal conduc- 
tivity measurements in the mixed state as a function of B 
give rise to additional information on low-lying excitation 
spectrum induced around vortices, which sensitively re- 
flects the gap structure In fact it has been already 
demonstrated that the orientational Sommerfeld coeffi- 
cient 7(0) and thermal conductivity k(9) show a charac- 
teristic oscillation with the relative angle 9 between the 
field direction and the direction of the gap node (or min- 
imum) H00HHII1E3- This allows us t0 P recisel y 
determine the gap anisotropy. 



As for observation of low-lying quasi-particle excita- 
tions in the mixed state, scanning tunneling microscopy 
(STM) provides a direct way to probe the spatial exten- 
sion of these quasi-particles STM experiments on 
2.ff-NbSe 2 comprise rich information on the gap struc- 
ture of this system [l4|. However in this method, as well 
as angle-resolved photoemission spectroscopy (ARPES) 
which yields momentum space information Hfil Il7j , 
observable materials are restricted due to the surface- 
sensitivity. 

Thus it is quite desirable to establish another method 
to probe the gap structure through bulk measurements. 
Here we propose the field-dependence of specific heat 
measurement at low temperature, namely the Sommer- 
feld coefficient j(B) in the mixed state. We demonstrate 
that the precise j(B) measurement at low field yields an 
important piece of information on the degree of the gap 
anisotropy, or the ratio of the minimum gap and maxi- 
mum gap in a superconductor. 

This type of experiments has been done in the past on 
a varietymaterials such as Nb 0, 2#-NbSe 2 [ml2c|. 
CeRu 2 0, LuNi 2 B 2 C 0, YBa 2 Cu 3 7 0, Sr 2 Ru0 4 
[l2| and Pr2- x Ge x CuOi-s 0. However, because of 
lack of theoretical detailed calculations, these data re- 
main largely to be analyzed. Theoretically j(B) was in- 
terpreted in terms of "rigid normal core concept" |25l |. 
that is, 7(B) oc B/B C 2 for superconductors with full 
gap because each flux carries a certain zero-energy den- 
sity of states (ZDOS). Volovik succeeded in describing 
qualitatively j(B) (cx y/B) for the nodal gap supercon- 
ductor [2fij |. However, there are many experimental data 
that j(B) shows convex curve, including y/B-\ike behav- 
ior, even in full-gap superconductors. So far no one has 
studied this problem in a quantitative level which allows 
us to examine abundant experimental data 7(B) in or- 
der to deduce the anisotropy of the gap structure, or the 
ratio of the minimum and maximum gaps, including the 
nodal case where the minimum gap is zero. 

The quasi-classical Eilenberger framework [23 ^ s mos t 
suited for the present purpose, which is microscopic valid 
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for almost all superconductors where ^> 1. We 
assume a superconductor in the clean limit and with 
the cylindrical Fermi surface, which does not alter our 
main conclusion essentially valid for three-dimensional 
case too. 

We introduce the pair potential A(r, 9), the vector 
potential A(r) and the quasi-classical Green's functions 
g(iu! n ,r,6), f (hj n ,T,9,) and p(iui n , r, 0), where r is the 
center of mass coordinate of the Cooper pair. The di- 
rection of the relative momentum of the Cooper pair, 
k = k/ |k|; k = (cos 9, sin 9), is denoted by the polar an- 
gle 9 relative to x-direction. The Eilenberger equation is 
given by 

L> n + ^v F • (j- + ^ A ( r )) | f(to>n,r, 9) 

= A(r,%(iu; n ,r,0), (1) 

{«n-^.(I-^A(r))}/t(W,r,tf) 

= A*(r,%(iw n ,r,0), (2) 
g(iu n ,r, 9) = [1- /t(L; n , r, 0)fQw n ,r, 9)]i, (3) 

where Reg(ko n , r, 9) > and vp = i»Fk is the Fermi ve- 
locity. <po is flux quantum. 

The external field H is applied along z-direction. With 
the symmetric gauge, the vector potential is written as 
A(r) = (1/2)H x r + a(r). Here a(r) is related to the 
internal field h(r), where h(r) = V x a(r). The pairing 
interaction is assumed separable V(9,9') — Vo<f>(9)(f>(9') 
so that the pair potential is A(r, 9) = A(r)</>(0). 

We numerically H H 13 solve Eqs. ©-© with the 
self-consistent conditions for A(r, 9) and a(r); 

A(r) = N V 2irT / f-<P(9')f^n,r,9'), (4) 
J W-^^ I ^^,r,9), (5) 

where j(r) = V x V x a(r), N is the density of states 
at the Fermi energy in the normal state. The cut-off 
energy is set as w c = 20T C . k = v /7C(3)/72(A /T c )k gl . 
C(3) is Riemann's zeta function. The Ginzburg-Landau 
parameter is set as kql = 9.0. Ao is the uniform gap at 
T = 0. 

The local DOS of the energy E is given by 

/■27T in 

N(E,r)=N j^ —Reg{iw n ^E + ir],r,0), (6) 

where g is calculated by the above Eilenberger equations 
with iu) n — > E + ir], using the solution A and A(r) by the 
self-consistent calculation. We typically use r] = O.OlAo. 

The calculation for Green's functions is performed at 
T/T c = 0.1 within the vortex lattice unit cell, which is 
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FIG. 1: (a) and (b) show field dependences of ZDOS. The 
filled (open) circles are for the total (extend) ZDOS. In (c) 
and (d) density maps of ZDOS are shown around a vortex 
core at the center and neighbor vortex cores. B/B C 2 = 0.24. 
The density height of the vortex center is truncated by 1. 
Isotropic gap case ((a) and (c)) and line- node gap case ((b) 
and (d)) are presented. 

divided into 81 x 81 mesh points. In this study, we as- 
sume that vortices form a triangular lattice. One of the 
nearest neighbor vortices is situated along x-direction. 

The gap function is taken as (f)(9) = 1 for the isotropic 
gap case, and \<j>(6)\ — VI ~ cos 69 = y2 |cos3#| for the 
line-node gap case. The gap function with finite mini- 
mum gap is chosen as 4>(9) = (1 — acos60)/i/l + a 2 /2, 
where a denotes the degree of the gap anisotropy, whose 
sign and magnitude characterize the gap structure. Here 
hexagonal crystal symmetry is assumed. 

The total ZDOS or N(B) is the spatial 
average of N(E — 0,r), which is given by 
N(B) — J N(E = 0, v)dr/ J dr. This is related to 
7(B) at T — > limit. To characterize delocalized 
contribution to total ZDOS, wc define "extended ZDOS" 
by N cxt (B) = J C N(E = 0,I)dl/ J c dl, where the path of 
the integration J ■ ■ ■ dl is along the boundary of the 
Wigner-Seitz cell. Thus N ext (B) is a useful quantity, 
related to the states outside of the vortex core. This is 
directly measurable by STM, and would somehow be 
related to the thermal conductivity. 

Let us now show our results deduced from extensive 
self-consistent calculations. We first take up two ex- 
treme cases for the isotropic gap and line node gap: In 
Fig. 1 (a), the field dependence of the ZDOS for the 
isotropic gap is shown. It is seen that N(B) is linear 
for smaller induction field B up to a certain field B*. 
Beyond B* N(B) becomes curving up to the upper crit- 
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FIG. 2: (color) (a) Field dependence of N(B) for anisotropic gap cases. Plotted data are for a = 0, 0.2, 0.3, 0.5 and the line- 
node from bottom to top. (b) The same data are replotted in logarithmic scale, (c) Field dependence of derivative 8N(B)/SB, 
i.e., slope of N(B). From right to left shadow lines indicate each B* for a — 0,0.2,0.3,0.5. In the inset, B* / B C 2 is plotted as 
a function of A m i n /A max . A m i n = <j)(0 = 0) (A max = 4>(9 = t/6)) is the minimum (maximum) of the gap function 4>(0). 



ical field B c2 . As for the extended ZDOS N eKt (B), it is 
concaved all the way up to B c2 . 

The presence of B* can be understood as follows: As B 
increases, considerable overlapping of ZDOS at each core 
makes the linear behavior in N(B) change to non-linear 
at B* . Up to B* , vortices carry independently a certain 
ZDOS to add up. We see it from Fig. 1 (c) that ZDOS is 
quite small in a large part of unit cell except for the small 
core region. This is compared with that for line node case 
in Fig. 1 (d). It is emphasized that the low field linear 
behavior in N(B) superficially accords with the so-called 
rigid normal core model |25l |. but the situation is more 
interesting and subtle: (1) There exists a crossover field 
B* from linear to non-linear even in the isotropic gap 
case, which is a key fact important in the anisotropic case 
below. (2) Even in small B region N ext (B) is substantial, 
meaning that vortex cores are overl app ed there. This 
invalidates rigid normal core picture [25j. 

The line node gap case as an extreme anisotropic 
one is shown in Figs. 1 (b) and (d). Both curves of 
N(B) and N cxt (B) are convex with no linear part at 
low B. These results are totally different from those 
of the isotropic case. In particular, the small field be- 
havior is affected by the presence of nodes. In the 
nodal gap case N ext (B) increases rapidly with B as seen 
from Fig. 1 (b), because the quasi-particle wave func- 
tions with zero-energy extend to far outside of the vortex 
core due to the node (see Fig. 1 (d)). Therefore N ext (B) 
shows the convex behavior at lowest B, which is con- 
trasted with the concave behavior of N cxt (B) at small 
B in Fig. 1 (a). Moreover N(B) is non-linear from the 
lowest field because the quasiparticles near the gap node 
are free from confinement in the core region, as clearly 
seen in Figs. 1 (b) and (d). 

These results again seem to coincide superficially with 



Volovik's Doppler shift calculation 26]. Here we em- 
phasize that Volovik calculation neglects the core con- 
tribution coming from the localized quasi-particles, tak- 
ing into account the extended quasi-particle contribution. 
Volovik result (N(B) oc ^/B) is, strictly speaking, only 
valid for smallest B region and is qualitative in nature. 

We now come to our main results for general 
anisotropic gap case; In Figs. 2 (a) and (b), we display 
N(B) for various values of a, including the above extreme 
cases with line node. It is seen from these that (A) For 
a ^ 1 N(B) is linear in B at lower fields. (B) The lin- 
ear B region becomes limited to lower B as a increases. 
This can be easily seen in logarithmic plot of Fig. 2 (b) , 
where the slopes (a ^ 1) in logarithmic plot change from 
N(B)/N„ cx B/B c2 to (B/B c2 ) 05 as B increases. There 
exists always a crossover field B* for general a, whose 
precise value is determined below. (C) At higher field 
near B c2 , N(B) behaves similarly, independent of a val- 
ues. It is clear that N(B) is not described by a single 
exponential form N(B) oc B p covering all B region ex- 
cept for the nodal gap case. Near B c2 the exponent p is 
approximated to be 0.5 for all curves. 

These results are interpreted along the same line above, 
namely, in terms of quasi-particle's localized-delocalized 
crossover upon varying B. In the region B < B* each 
vortex core independently contributes to the DOS at zero 
energy while B > B* the quasi-particles are brought 
with flux lines overlapping each other, making N(B) non- 
linear. Thus a question is how B* is correlated to the gap 
anisotropy a. 

We exhibit the derivative SN(B)/5B in Fig. 2 (c). If 
N(B) is linear in B, the derivative is a constant and 
equals to the slope of N(B) near B = 0. This linear slope 
constant, indicating how much DOS is brought into a 
system by a flux line, increases with a. This implies that 



4 







1 


■ 


• 


2//-NbSe 2 






(Hanaguri) 




o 


Full gap (a = 0.5) 


I i 


A 


Line-node gap 



B/B C 2 

FIG. 3: (color) Comparison with the 7(B) in 2#-NbSe 2 
with our results for a — 0.5 and line node gap. 



exists the Fermi velocity anisotropy, one may concern its 
effect on the present conclusion. It turns out that it does 
not alter our conclusion essentially [13 ■ (4) Although we 
considered here 2D system, 3D effect also is checked. We 
find that it does not also change our conclusion because 
N(B) is a bulk quantity. Thus we believe that a care- 
ful N(B) measurement at low enough T is a powerful 
method to identify the gap anisotropy. 

We are grateful to Y. Matsuda, T. Sakakibara, K. 
Izawa, T. Hanaguri, T. Kiss and J. E. Sonier for use- 
ful discussions. One of the authors (N. N.) is supported 
by a Grant-in- Aid for the 21st Century COE "Center for 
Diversity and University in Physics" . 



the ZDOS per vortex increases with the gap anisotropy. 
Physically it is because the average core radius measured 
by the extension of ZDOS around a core is effectively 
larger for larger gap anisotropy. 

It is clearly seen from Fig. 2 (c) that there is a con- 
stant 5N(B)/SB region up tp B* for each value of a. 
As shown in the inset of Fig. 2 (c), this crossover field 
B* progressively decreases as the ratio of the minimum 
and maximum gaps grows, or a increases towards 1. In 
the nodal gap case, the derivative 5N(B)/SB is diverg- 
ing because the exponent p in N(B) is less than unity 
[28j. Thus it is possible to estimate the degree of the 
gap anisotropy a by identifying B* through specific heat 
measurement. 

For illustrative purpose, we show the comparison be- 
tween the 7(B) data of 2#-NbSe 2 and our results 
(a = 0.5 and line node case) in Fig. 3. 2iJ-NbSe 2 is a 
anisotropic but full gap superconductor. While it is pos- 
sible even to say that 7(B) cx \fB at high fields, j(B) de- 
viates from behavior at low fields. The low field data 
fit to a = 0.5 rather than line node gap. The anisotr opy 
a = 0.5 is also estimated by the recent ARPES data 17] 
and also by our previous identification |29| . 

In order to identify the a value in the actual measure- 
ment, we must examine the applicability of the present 
approach. The possible obstacles, which are not con- 
sidered here are following: (1) Non-magnetic impurities 
certainly modify N(B) behavior as discussed elsewhere 
[3fl l3ll ] . One needs samples in clean limit. (2) The actual 
superconductors are often multi-band system, or multi- 
gap one. Typical example is MgB 2 where distinct two 
different gaps are assigned for two-bands. In fact N(B) 
consists of two linear fines J32J, ea ch segment is explained 
in terms of isotropic gap |33l l34| . Thus it is possible to 
separate out the multi-gap effect. Such multi-gap be- 
havior is seen in various systems, e.g. Sr2Ru04 [l2j |. 
2£/-NbSe 2 H3 or CeRu 2 EJ etc. For 2iJ-NbSe 2 , other 
measurements [lH imply multi gaps. (3) Since there 
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